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We consider spherically symmetric distributions of anisotropic fluids with a central vacuum cavity, 
evolving under the condition of vanishing expansion scalar. Some analytical solutions are found 
satisfying Darmois junction conditions on both delimiting boundary surfaces, while some others 
require the presence of thin shells on either (or both) boundary surfaces. The solutions here obtained 
model the evolution of the vacuum cavity and the surrounding fluid distribution, emerging after a 
central explosion, thereby showing the potential of expansion-free condition for the study of that 
kind of problems. This study complements a previously published work where modeling of the 
evolution of such kind of systems was achieved through a different kinematical condition. 
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I. INTRODUCTION 

Highly energetic explosions in self-gravitating fluid dis- 
tributions are common events in relativistic astrophysics 
(see for example p], [|J and references therein) . There are 
two main questions related to this issue: 

• What is the energy source of such explosions? 

• How the system evolves after the explosion? 

This work is devoted to the second of the above ques- 
tions. Our approach is completely analytical, which 
presents both advantages and disadvantages. On the one 
hand the analytical approach allows one to carry out a 
rather model independent analysis, allowing in turn to 
obtain general conclusions about the pattern evolution 
of the system. On the other hand however, in order to 
handle analytically the involved expressions, one needs to 
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introduce simplifying assumptions which might exclude 
important physical phenomena. 

We are interested in analytical models which even if 
might be relatively simple to analyze, still contain some 
of the essential features of a realistic situation. It should 
be emphasized that we are not interested in the dynamics 
and the conditions of the creation of the cavity itself, but 
only in its evolution once it is already formed. 

A similar study to the one presented here was recently 
published [3j, however the main difference between the 
two being the kinematical condition imposed on the evo- 
lution of the system. In Q it was assumed that the 
proper radial distance between any two infinitesimally 
close fluid elements remains constant all along the evo- 
lution (the purely areal evolution condition). Here we 
assume that the expansion scalar vanishes for all fluid el- 
ements. As discussed in either assumption creates 
conditions for the existence of a vacuum cavity surround- 
ing the centre. The dynamical instability of expansion- 
free fluids is studied in @. 

We have two hypersurfaces delimiting the fluid. The 
external one separating the fluid distribution from a 
Schwarzschild spacetime (for simplicity we assume the 
evolution to be adiabatic) and the internal one, delim- 
iting the cavity within which we have Minkowski space- 
time. It should be clear that for cavities with sizes of the 
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order of 20 Mpc or smaller, the assumption of a spheri- 
cally symmetric spacetime outside the cavity is quite rea- 
sonable, since the observed universe cannot be considered 
homogeneous on scales less than 150-300 Mpc. Never- 
theless for larger cavities it should be more appropriate 
to consider their embedding in an expanding Lemaitre- 
Friedmann-Robertson- Walker spacetime (for the specific 
case of void modeling in expanding universes see 0, H| 
and references therein). 

Thus, we have to consider junction conditions on both 
hypersurfaces. Depending on whether we impose Dar- 
mois conditions Q or allow for the existence of thin shells 
10], different kinds of models are obtained. In this paper 
we are particularly interested in models satisfying Dar- 
mois conditions, however some models presenting thin 
shells will be also exhibited. 

The first known solution satisfying the vanishing ex- 
pansion scalar is due to Skripkin [ill, EH- ^ is worth 
noticing that Skripkin did not assume explicitly the ex- 
pansion scalar to vanish, however the conditions he im- 
posed (isotropic fluid with a constant energy-density) 
imply the vanishing of the expansion scalar. Unfortu- 
nately though, as it was shown in [5] , the Skripkin model 
does not satisfy Darmois junction conditions (in principle 
Skripkin models with thin shells satisfying Israel condi- 
tions can be obtained). 

It is our purpose here to find exact solutions under 
the vanishing expansion scalar condition and satisfying 
Darmois junction conditions on both delimiting hyper- 
surfaces. 

For sake of generality we shall consider an anisotropic 
fluid. Indeed, as it has been shown in [4], in the case of 
isotropic pressure, the expansion-free conditions imposes 
that the energy-density is independent on the time-like 
coordinate, which severely restricts the models (further 
arguments to justify such kind of fluid distributions may 
found in |13l - [l5| and references therein). A detailed de- 
scription of this kind of distribution, the Darmois junc- 
tion conditions on both, the inner and the outer bound- 
ary surface, as well as definitions of kinematical variables, 
and basic equations, are given in section II. The general 
description of the expansion-free fluids is given in section 
III, whereas explicit solutions are found and described in 
sections IV and V. Finally results are summarized in the 
last section. 



II. FLUID DISTRIBUTION, KINEMATICAL 
VARIABLES AND BASIC EQUATIONS 

We consider a spherically symmetric distribution of 
collapsing fluid, bounded by a spherical surface Y,( e K 
The fluid is assumed to be locally anisotropic (principal 
stresses unequal). 

Choosing comoving coordinates inside £( e ) , the general 
interior metric can be written 



where A, B and R are functions of t and r and are as- 
sumed positive. We number the coordinates x° = t, 
x 1 = r, x 2 — 9 and x 3 = (p. Observe that A and B 
are dimensionless, whereas R has the same dimension as 
r. 

The matter energy-momentum T~p inside T,^ has the 
form 

T c7/3 = (M + P±)V a Vp + P^g a p + Tl Xa Xp, (2) 

where /i is the energy density, H = P r — P± , P r the radial 
pressure, P± the tangential pressure, V a the four velocity 
of the fluid, and x" a unit four vector along the radial 
direction. These quantities satisfy 

V a V a = -1, X a X a = 1, X°V a = 0, (3) 

since we assumed the metric ([1]) comoving then 

V a = A" 1 ^, x a = B^S?. (4) 



A. Einstein equations 

From (p} and ©, Einstein equations 

G a/ 3 — 8irT al 3 , 

read 

8nT 00 = SnfiA 2 = I 2- t 
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where the prime stands for r differentiation and the dot 
stands for differentiation with respect to t. 



B. Kinematical variables and the mass function 



ds 2 



-A 2 dt 2 + B 2 dr 2 + R 2 ( 



0, (1) 



The four acceleration a a is given by 

a a = V a . l0 V , 



(10) 
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from (ITU)) and ((H) we obtain 

■"■ 2 n 

A 



AB J 



with a a = ax a ■ The expansion is given by 



J ~ V " a ~ A \ B R 



and for the shear we have 



= V( Q;( 9) + a,( a Vp) - -Qh a p. 



(11) 
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where h a p = g a /3 + V a Vp. Using (|4| we obtain the non- 
vanishing components of (|13[) 
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cr Q « may also be written as 
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Next, the mass function m(t, r) introduced by Misner 
and Sharp [lj| (see also [l?}) is given by 
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We can define the velocity U of the collapsing fluid as 
the variation of the areal radius R with respect to proper 
time, i.e. 



Then (fT8l) can be rewritten as 
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With the above we can express ([7]) as 

0= J-(e-cr)'-^. 
3R' V 7 i? 

While from (TT5| we have 



m = -4wP r RR , 



and 



to' = AirpR'R 2 . 
Combining (|2"21 and (|2"3")) we obtain 

fiR' + P' r R + {P r + n) [ Rl + 2R 
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The integration of (|2"B")) produces 

4TrnR 2 R'dr 
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(22) 
(23) 

(24) 
(25) 



where a regular centre to the distribution, m(0) = is 
assumed. We may partially integrate (|25|) to obtain 
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3^=47r/i-^ / p'R 3 dr. 
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C. Junction conditions 

Outside £' e ) (whose equation is r = r e = constant) we 
assume the Schwarzschild spacetime, described by 

ds 2 = - f 1 - 2 — ) rfy2 - 2d /^ + P 2 ^ 2 + sin2 6,d< / ,2 )> 

(27) 

where M denotes the total mass and v is the retarded 
time. 

The matching of the adiabatic fluid sphere to 
Schwarzschild spacetime, on the surface r = r e = con- 
stant (or p = p(v) e in the coordinates of (|27[) ). requires 
the continuity of the first and second differential forms 
(Darmois conditions), implying (see [l8[ for details) 



(19) 



(20) and 



Adt = dv fl - 2— J , 

i? S = p{ V ), 
m(t,r) E =' M, 
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(30) 
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where = means that both sides of the equation are eval- 
uated on £( e ) (observe a misprint in equation (40) in [l8T ] 
and a slight difference in notation). 

Comparing (l3Tj) with ([7]) and (JHJ) one obtains 



P r = 0. 



(32) 



As we mentioned in the introduction, the expansion- 
free models present an internal vacuum cavity. If we 
call E'*) the boundary surface between the cavity and 
the fluid, described by the equation r = = constant, 
then the matching of the Minkowski spacctimc within the 
cavity to the fluid distribution, implies 



s (<) 

m(t,r) = 0, 

y(i) 
P r = 0. 



D. Weyl tensor 



(33) 
(34) 



Finally, comparing (|26|) with (j40|) we get 

47T f r 

£ = — fi'R 3 dr - 47rII, 
" Jo 



(41) 



which exhibits the dependence of the Weyl tensor on 
energy-density inhomogeneity and local anisotropy of 
pressure. 



E. Bianchi identities and an evolution equation for 
the Weyl tensor 

The two independent components of Bianchi identities 
for the system under consideration read (see [3| for de- 
tails): 



fi+ (fi + P r )AG- 2n^ = 0, (42) 
R 



The Weyl tensor is defined through the Riemann tensor 
^a/3u> t ne Ricci tensor R a p and the curvature scalar TZ, 
as 



+1^9*0 + \n{5^g ali - g aP SP). (35) 



The electric part of Weyl tensor is defined by (the mag- 
netic part vanishes due to the spherical symmetry) 



E a /3 = Ca^fivV^V ', 

with the following non- vanishing components 
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where 



(36) 



(37) 



P' r + ^ + P r )^ + 2H-| = 0. (43) 

Finally, the following evolution equation for the Weyl 
tensor will be used, which may be derived from the 
Bianchi identities [l9|, [2(| (see also [l4| for details) 

[£ - 4x(p - IT)]' + 3[£- 47r(/x + P ± )} J = 0. (44) 



III. THE EXPANSION-FREE FLUID 

If we assume the evolution to be expansion-free, i.e. 
O = 0, then from (fT2l) wc have 



B R h 

and after substituting (|4"5j) into 0, we get 
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Observe that we may also write E a p as 



E. _ — 

~R~ ~A 



where r(t) and h(r) are arbitrary functions of their ar- 
(38) guments. Thus, without loss of generality, we shall take 
h — 1 (observe that with this choice a unit constant with 
1 dimensions [r ] is assumed to multiply dr 2 ). 

Jjp ■ Next, using the equations ([22]l. pB]l. (|42|) and (|4"5]l. the 

physical variables /x, P r , and II, can be written through 
the mass function m and the metric function R as follows, 



E a p = £ ^XaXfi - \ h <*p \ ■ (39) 
Using ©, ©, © with (HD and ([35) we obtain 



3^5 = 4tt( M - H) - f . (40) 
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(47) 
(48) 
(49) 
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On the other hand, combining (TT5|) . (|45l) and f4"6")l we 

obtain, 



m(t, r) 



- - 

2 \R 3 



R' 2 R 5 + R 



(50) 



To completely determine the metric we need to give an 
extra condition. In the next sections we provide some 
options. 



A. a=0 



For a = (j = 1), integration of ([58|)produces. 



46 1 
3c ~ 8 



c{r + (i) 2 



(59) 



implying 



IV. SOLUTION I 

The first family of solutions will be obtained by choos- 
ing the function m(t, r) as follows (2lj . 

2m(t,r) = jR+-kR 3 + -IR 5 , (51) 
3 5 
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where /3(t) is an arbitrary function oft. Then, using (|52[) . 
(|48| and (|49|) the physical variables, fi, P r and IT, can be 
expressed as 



where j(t), k(t) and l(t) are so far arbitrary functions of 
t. Substituting (f5T|) into (|4"T)) we obtain for the energy 
density 
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87r / i= 777 +k + lR\ 
R z 



(61) 



87r/i= 4t +k + lR' 2 . 
R z 



(52) 



Then, taking the r derivative of ([5H)) and using (|4T1) and 
(l52l we obtain, 



1 + kR 2 + /i? 4 i? / k I 2 \ 
SnP r = - 2 ^ » , (62) 



j - 1 + fci? 2 + ^i? 4 + 3— + 2i?"i? 5 + 5i?' 2 i? 4 = 0. (53) 
Combining (|50]l. ([5T]) and ([53]). we get 
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A: + Zi? 2 
2RR 



(63) 



4(j - 1) + 2fci? 2 + -IR 4 + 2R"R 5 + 8R' 2 R 4 = 0, (54) 
5 

which by introducing the new function Z = R 2 , can be 
transformed into 

a + bZ + cZ 2 + Z"Z 2 + ^Z' 2 Z = 0, (55) 

where 



This subfamily of solutions may satisfy Darmois con- 
ditions. Indeed junctions conditions (|3U)) . (|3"2"j) and (|3"3"1) . 
(|3"4l lead to two independent equations for the three func- 
tions b(t), c(t) and (3(t) which can be satisfied by any 
convenient choice of one of these functions. However for 
all the attempts we have carried out resulting expressions 
were extremely cumbersome and therefore we shall not 
pursue this direction further. 



a(t)=4(j-l), b(t) = 2k, c(t) = -l. (56) 

5 

Finally, feeding back Y(Z) = Z' 2 into ([gB]l. leads to the 
first order equation 



B. b=c=0 



dY Y 

dz +3 z 



Z 2 Z 



and after integration we obtain 



(57) 



(58) 



For b(t) = c(t) — 0, integration of (|58|) produces 



2/3 



(64) 



or 



This last equation can be integrated in its general form, 
however in order to obtain simple models, in what follows 
we shall consider specific cases. 



R 



-ar + f3 



1/3 



(65) 
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Then, using ([52")). ((38"]) and (gHl) the physical variables //, 
P r and II read, 



-2/3 

*-/' = ./'( i^or + iS) ,(66) 



87rP r = 2i( -V=^r + /3 



1/3 



2 \ —a 3 



87rn=-i( -^r + /3 



/3 \ - 2 /3 

-jf-^r + /3) ,(67) 

1/3 /I d 2-^ 1 

; r P 

2 -J— a 3 



-2/3 

^r + /3) .(68) 



It is a simple matter to check that this subfamily of 
solution does not satisfy Darmois conditions on the in- 
ternal boundary surface r — r,. Accordingly a thin shell 
should be present there. 



V. SOLUTION II 

The next family of solutions is found with the addi- 
tional assumption P r — 0. Fluid spheres with tangential 
stresses alone, have a long and a venerable history start- 
ing with Lemaitre seminal paper [22| . Afterwards, many 
authors have used such kind of fluid distribution for dif- 
ferent purposes (see for example [23T - I291 ] and references 
therein). 

Thus, assuming P r = 0, we find from ((2"4"]) 



V = 



P'P 2 



=> R = 3 



dr + C 2 , 



(69) 



where C\ (r) and Ci (t) are two functions of integration. 
Using (|47p in (]69p . the former can be easily identified as 



_ m 
6l " 4^ 



(70) 



Now, feeding back (l4"5j) in (|4"2")) we obtain for the tan- 
gential pressure P± 



P 



fi R 
2R' 



(71) 



In what follows we shall impose different additional re- 
strictions in order to obtain different subfamilies of mod- 
els. 



Assuming 



P r = 0, Pi 



afi, 



(72) 



where a is a constant, then (|4"3"]) becomes 

R' „ / , \ R' 
— = 2(a- 1 — , 
R K 1 R' 

and after integration 

R = fR 2 ^ a - 1 \ R'=gR 2 ^ 



(73) 



(74) 



with f(t) and g(r) denoting arbitrary functions of their 
arguments. Then from (|74l) we obtain, 



P 



3-2a 



(75) 



with 



m = (3 - 2a) J fdt, X (r) - (3 - 2a) J gdr. (76) 

Without loss of generality we may choose r(t) — f, im- 
plying because of (|4"6")l and (1741 . 

A = R 2a . (77) 
Considering (|72|) with (|22f) we have 

m = m(r), (78) 
and from ([18]) with (J45J) and <j77j) we obtain 
P / P 2 



m{r) = f I P^ " ^ + 1 



(79) 



Matching conditions (|3U|) and (f3"3"() on the hypersurfaccs 
r = r e and r = r,-, together with ((79"]) . produce 



and 



P 2S = R 4a (g 2 R 4a - 1), 



£2 SW ^ 4q f ^£ , „2 d 4q 



,g z P 4 " - 1 



(80) 



(81) 



Equation ([8T)]) can be integrated for any given value of 
a (the particular case of dust, a — has already been 
considered in [5j). 

For a = 1/2 the solution reads 



P 



3 W 1 



gcos(t + to) ' 
and from ((75)1 evaluated on £W, 



1 



<? 2 cos 2 (i + < ) 



(82) 



(83) 



Thus the time dependence of all variables is fully deter- 
mined. The radial dependence (Ci or x) c & n be obtained 
from the initial data. 
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However, this solution is inconsistent for any a > 
with junction conditions satisfied on both hypersurfaces. 
Indeed, from ([50)1 and (|5T1) evaluated at t = we get 



Now substracting (|85p from (|84p we obtain 



and 



i?r(o) 
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I+ I + 4. 9 ^(0) 
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i?e(0) 



Re(0) 



Rt a (0) 



25 2 



1 -J^r] + 4 5 2 i? 2 (0) 



1/2' 



>(85) 



i?^(0)-i?f(0) = — 



2 i? 2 (0) 



1/21 
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I - 



2M 

i?e(0) 



I - 



2Af ~ 



1 1/2' 



+ - 



2 i? 2 (0) 



(86) 



which should be negative for all t (if a > 0) (or zero if 
a = 0) since obviously 72 e (t) > Ri(t). Unfortunately this 
is not the case. 

In fact, from the definition of U, using (|4"6]) and ([77)1 
we obtain 



1 R] Rl a ' 



(87) 



and 



Rl Rl« 

Since a < 1 and R e > Ri it follows from the above that 

Ri > R e . (89) 

Then using {89]) in ([86]) we obtain at once that Rf a (0) > 
i?g Q (0) which is unacceptable. 

Therefore such models require the presence of thin 
shells at either (or both) boundary surfaces. 



_ 3 fJ.op.oCi ( r 3 + ii C 2 
± ~ 2 r 2 \ ,v^3 ,,2/ 



,/i r 3 - [j^C-2 
and from (|50|) and (f9"Tj) we obtain 



i 



(92) 



(93) 



Evaluating this expression on £W by using (f3"3"]) we obtain 

„4 x -i 



J? 



4_ r 2 



Mo 



1 



(94) 



implying the restriction r\ > for all t. 

Next, from the definition of U and (|46|) we get 



U 



R?' 



(95) 



implying 



B. 



/i = fJoCi/r 



Next, the following subfamily of solutions will be ob- 
tained under the assumption that the energy-density is 
separable. Thus we write 



MoCi 



(90) 



where no(t) is an arbitrary function of its argument and 
d(r) is given by ([70]). Then substituting flSC) into (JBSJ) 
and ([7T]) . we get 



7? 



— + C 2 
Mo 



1/3 



(91) 



U 



Mo 



(96) 



Thus, from the absence of superluminal velocities (U < 
1) and ([TO]) we need to impose 



V2fj > rl > hq. 
Hence from ([95]) we have 



(97) 



(98) 



and since R e > Ri then from (f98|) it follows that the inner 
boundary surface moves faster than the outer one. 
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Now, since r is a function of t which may be chosen 
arbitrarily, we put without loss of generality 



implying 



At 



i Ri = Ui, T = Ria t = U e R 2 e 



(99) 



(100) 



Let us now consider the junction conditions at the 
outer boundary surface Using USD), dSSJ) and ([TOO]) 

we obtain 



(101) 



R e = 2M [U- - — + 1 
imposing the restriction on r e , 



(102) 



which is the same as imposed on r,. 

Thus if we provide the time dependence of /i then (|94[) 
allows to determine the time dependence of Ri , and using 
this result in (|9T|) we obtain C^it). Finally evaluating 
flgTJ on we get R e . 

Alternatively we may provide the time dependence of, 
either Ri or R e which allows also to find the time depen- 
dence of all variables. We shall next exhibit a specific 
model obtained from the additional restriction of confor- 
mal flatness. 



1. Conformally flat models 

Instead of assuming any specific time dependence for 
any of the variables considered above, we shall assume 
here that the fluid is conformally flat, i.e. we shall assume 
8 = then we obtain from (J40J) (with P r = 0) 



H + P± = 



3m 



(103) 



Next, evaluating (fTU3")l at EW, using $TQ) and 0, we 
obtain 



(104) 



Ri = a/iQ, 



where a is a constant. Feeding back (|104[) into (|96|) we 
obtain 



(105) 



with fi 2 = rfaf. Thus 



R,. = ± 



1/2 



(106) 



where + (— ) corresponds to expanding (contracting) con- 
figurations. The integral of (|106|) is expressed through 



elliptic functions, and we shall not display explicitly the 
solution. 

The time dependence of the remaining functions is ob- 
tained as in the previous case. 

Finally, we need to determine C\ (r) , for doing that we 
may proceed as follows: from (|2U)) and (|M|) we have 



47T 



R' R 2 HqCi 



(107) 



Taking the r derivative of m from (|93|) . feeding it back 
into (|107[) and evaluating the result at t — we obtain 
Ci(r). 

These solutions are regular everywhere (within the 
fluid distribution) and satisfy basic physical restrictions 
such as positivity of energy density distribution and 
M>-P±. 



VI. CONCLUSIONS 

We have revisited the original Skripkin problem, but 
now allowing for a more general fluid distribution. This 
enabled us to find analytical models which satisfy Dar- 
mois junction conditions on both delimiting boundary 
surfaces. Obviously, relaxing Darmois conditions on ei- 
ther (or both) boundary hypesurfaces enlarges the fam- 
ily of possible solutions, some of them are explicitly ex- 
hibited in previous sections (for models of voids within 
the thin wall approximation, see [30|-(32l| and references 
therein). 

The most important contribution in this manuscript 
is to illustrate the advantages of using the expansion- 
free condition for the modeling of cavity evolution. For 
that purpose we have integrated analytically the perti- 
nent equations under a variety of circumstances, proving 
that such an integration could be achieved without un- 
due difficulty. The obtained solutions were not intended 
to describe any specific astrophysical scenario but just to 
bring out the potential of the expansion-free condition to 
model situation where cavities are expected to appear. 

Possibly our solutions could be applied as toy models 
of localized systems such as supernova explosion mod- 
els. Most likely the combination of the vanishing ex- 
pansion scalar condition with numerical integration of 
corresponding equations would produce more physically 
meaningful models. 

It is worth mentioning the fact that in spite of the van- 
ishing of the expansion scalar, the energy-density changes 
with time. This effect is due to the work done by the 
anisotropic force as it appears from (l4"2j) . 



Acknowledgments 

LH wishes to thank Fundacion Empresas Polar for fi- 
nancial support. ADP acknowledges hospitality of the 
Departamento de Fisica Teorica e Historia de la Fi'sica 



9 



at Universidad del Pais Vasco. JO wishes to thank Uni- No. FK10. 
versidad de Salamanca for financial support under grant 



[1] A. Mazure and S. Basa Exploding Superstars, (Springer, 
Berlin) (2009). 

[2] G. Ghisellini arXiv:astro-ph/1010.3015vl (2010). 

[3] L. Herrera, G. Le Denmat and N. O. Santos Class. Quan- 
tum Grav. 27, 135017 (2010). 

[4] L. Herrera, N. O. Santos and A. Wang Phys. Rev. D 78, 
084026-10, (2008). 

[5] L. Herrera, G. Le Denmat and N.O. Santos Phys. Rev. 
D 79 , 087505, (2009). 

[6] L. Herrera, G. Le Denmat and N. O. Santos arXiv.gr- 
qc/1103.3617 (2010). 

[7] W. Bonnor and A. Chamorro Astrophys. J. 361, 21 
(1990). 

[8] R. Torres Class. Quantum Grav. 22, 4335 (2005). 

[9] G. Darmois Memorial des Sciences Mathematiques 

(Gauthier-Villars, Paris, 1927) Fasc. 25. 
[10] W. Israel II Nuovo Cimento B 56, 1 (1966). 
[11] V. A. Skripkin Soviet Physics-Doklady 135, 1183 (1960). 
[12] H. Stephany, D. Kramer, M. MacCallum, C. Honse- 

laers and E. Herlt Exact Solutions to Einstein's Field 

Equations. Second Edition, (Cambridge University Press, 

Cambridge) (2003). 
[13] L. Herrera and N. O. Santos Phys. Rep. 286, 53 (1997). 
[14] L. Herrera, A. Di Prisco, J. Martin, J. Ospino, N. O. 

Santos and O. Troconis Phys. Rev. D 69 084026 (2004). 
[15] B. Ivanov, Int. J. Theor. Phys. 49, 1236 (2010). 



[16] C. Misner and D. Sharp, Phys. Rev. 136, B571 (1964). 
[17] M. Cahill and G. McVittie, J. Math. Phys. 11, 1382 
(1970). 

[18] R. Chan Mon. Not. R. Astron. Soc. 316, 588 (2000). 
[19] G. F. R. Ellis, Relativistic Cosmology in: Proceedings 

of the International School of Physics " Enrico Fermi", 

Course 47: General Relativity and Cosmology. Ed. R. K. 

Sachs (Academic Press, New York and London) (1971). 
[20] G. F. R. Ellis, Gen. Rel. Grav. 41, 581 (2009). 
[21] M. Chaisi and S. D. Maharaj Gen. Rel. Grav. 37, 1177 

(2005). 

[22] G. Lemaitre Ann. Soc. Set. Bruxelles A 53, 51 (1933). 
[23] A. Einstein Ann. Math. 40, 4 (1939). 
[24] B. K. Datta Gen. Rel. Grav. 1. 19 (1970). 
[25] H. Bondi Gen. Rel. Grav. 2 321 (1971). 
[26] P. S. Florides Proc. R. Soc. London A 337, 529 (1974). 
[27] L. Herrera and N. O. Santos Gen. Rel. Grav. 27, 1071 
(1997). 

[28] G. Magli Class. Quantum Grav. 15, 3215 (1998). 

[29] S. M. C. V. Gongalves, S. Jhingan and G. Magli Phys. 

Rev. D 65, 064011 (2002). 
[30] K. Lake and R. Pim Astrophys. J. 298, 439 (1985). 
[31] R. Pim and K. Lake Astrophys. J. 304, 75 (1986). 
[32] R. Pim and K. Lake Astrophys. J. 330, 625 (1988). 



